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Abstract. 

The dynamical eflFects on the scale factors due to the scalar (/)-field at the early stages 
of a supposedly anisotropic Universe expansion in the scalar-tensor cosmology of 
Jordan-Brans and Dicke is studied. This universe shows an isotropic evolution but, 
depending on the value of the theorie's coupling parameter u), it can begin from 
a singularity if u; > and after expanding shrink to another one; or, ii uj < 
and —3/2 < uj < —4/3, it can evolve from a flat spatially-inflnite state to a non 
extended singularity; or, if —4/3 < u; < 0, evolve from an extended singularity to a 
non singular state and, at last, proceed towards a singularity. 
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Nowadays, and despite the lack of favorable observational evidence in todays uni- 
verse, there is a great deal of interest in the Jordan-Brans-Dicke scalar-tensor theory 
of gravity (JBD) [1, 2] — and in other scalar tensor theories too — because of the 
emergence of superstring theories which lead naturally to a dilaton theory of grav- 
ity where scalar fields are mandatory, and also due to the emergence of extended 
inflation models and pre-big bang cosmologies where scalar fields can provide solu- 
tions to some of the problems of inflation [3-5]. This comes about since the JBD 
action functional already includes a string sector where the dilaton fleld (f) d can be 
suitably related with the JBD scalar fleld cf) oc exp(— ^i:)). The important role of 
the 0-fleld of JBD would especially occur at the strongly relativistic stages of the 
Universe [5, 6]. Thence, the importance of studying the early phases of the homo- 
geneous and anisotropic Bianchi universes in JBD cosmology, since is possible that 
the early cosmological expansion of the universe be determined by an anisotropic 
but homogeneous vacuum stage with a non vanishing scalar fleld. 

In this work we report the early stages of evolution of the Bianchi type VIIq 
vacuum universe and show that, despite the suppposed anisotropic behaviour of 
the Bianchi universes, it evolves isotropically. The beginnings of such universe 
are shown to depend on the value of the coupling parameter uj appearing in the 
JBD theory. If the coupling parameter is positive u > 0, the universe starts from 
a singularity, expands reaching a maximum size and then contracts until another 
singularity is reached. If the coupling parameter is instead negative, then the 
behaviour depends on whether u is larger or smaller than the critical value Uc — 
—4/3. In the former case (0 > a; > Uc), the universe begins from a spatially-inflnite 
singularity, it quickly collapses to a flnite state with non zero curvature and then, 
more slowly than in the previous stage, reaches again a singularity. In the latter 
case {uJc > UJ > —3/2), on the other hand, the universe starts in a spatially inflnite 
but flat state and, from there, it collapses to a non-extended singularity. 

To obtain the Bianchi- type VII equations, let us write the line element of the 
spacetime, using signature +2 and natural units c — G — 1, as 

ds^ ^ -dt^ + hij{t)u!'u\ (1) 
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where the hij (t) is the 3-metric on the surface of homogeinity, t is the synchronous 
or cosmological time, a;* are the one-forms [7]: 

W"*" —ai {{i] — kv)dy — vdz) , 
(J^ —a2 {fdy — [t] -\- kv)dz) , 
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(jj —a^dx, 
^dt, 

where rj — exp(— /ca;) cos(Ma;), i/ — {—M~^)eKp{—kx)sm{Mx), k — h/2 and 
M — {1 — k'^Y^'^ ; the parameter h distinguish the case analysed here {h — 0) from 
the generic singular with h 0. If we now insert the line element (1), with the 
forms (2), into the JBD field equations in a vacuum [2, 8, 9], we get 



|^(lnaO + I (In a.) I (In aia^aa) + |(lnaO i ^ ^ ^^^^^ 



+ + ^ip2 + :Fip3 = 0, 

where one of the coordinates has been chosen as the synchronous time t. We 
additionally have what we have called the constriction equation, coming from off- 
diagonal terms in the JBD field equations. 



|(lna0|(lna,) + |(lnai)|(lna3) + |(lna2)|(lna3) 

/ • \ / ; \ ^ (4) 
+ ^(Inaiasaa) f |) " | (| ) + -^^""i^ + ^^^^ + ^^^^ = 0' 

finally, the JBD scalar field comply with 



|([aia.a3]^)=0, (5) 



where /3i = [ai/ {2ajak)Y ^^'^ the indexes z,j,/c are to be taken in cyclic order of 
1,2,3. Equations (3), (4) and (5) are written in the standard form we introduced 
previously for solving the Bianchi models in JBD [10]; the specific values for the 
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constants appearing in them are: Ai — 4M^ — (5/2), A2 — Ai — 2, As — 0, 
A4: — Ai — 1, Si — S3 = J^i — J^2 — —2, S2 = J-'s — 2, S4 — — —1, and 
= 1 — (/i^/4). Notice that according to these relationships, h must be restricted 
to be \h\ < 2. The equations for the Bianchi VIIq model can be obtained from 
these just by taking h — 0. The Bianchi models, as consequence of certain features 
of their metrics, lead to additional relationships between their scale factors. In this 
case, we got two additional relationships [9] 

|(lna,)-^(lna.) = 0, (6) 

and 

= 0. (7) 

Equation (6) always provides an additional relationship between the two scale fac- 
tors, ai and 02, not mattering what the value of the /i-parameter; on the other 
hand, when /i 7^ 0, equation (7) leads to a singular universe in which the surfaces of 
homogeneity collapse to a 2-manifold, to a one-dimensional manifold or, even, to a 
single point. A more detailed analysis of this case is reported in [9, 11]. Equation 
(7) becomes just a trivial identity when h — 0, not restricting in any way the values 
of the scale factors. This is precisely the case analysed in this contribution. 
Prom equation (5) we easily obtain 

d(f) 

ai ^2^3 = 90, (8) 

where 0o is an integration constant. Thus, we can define the scaled 0-field, as 
$ = 0/00, also called the intrinsic time [9, 10], and we get the relationship 



dt = (010203) (9) 

We use $ as a new time coordinate in susbstitution of the cosmological time t for 
solving equations (3); $ has been found useful for analysing the Bianchi vacuum 
models in several situations [9, 10, 12, 13]. Let us introduce the notation ( )' = 
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and, if we define the Hubble expansion rates as = (Inaj)', the field equations 
become 

Hi+ ^ + Jiai + ]Ci4 + Afialal ^ 0, 1,2,3 (10) 

and the constriction equation becomes 

(Inaia2a3)' oo 



Hi H2 + Hi H^ + H2 H^ + 



$ 2$2 (11) 

+ J'^a^ + /C4 03 + A/4 02 03 = 0, 
the specific values for the constants appearing in equations (10) and (11) are 
combinations of the constants previously used: J^i — J^3 — ]Ci — ]C2 — —1/2, 
= /C3 = 1/2, J4 = /C4 = -1/4, M = 4M2 - (5/2), Ml- 2, A/'3 = 0, 

A/4 = A/i — 1. We have written the equations as to emphasize the relationship with 
our previous work on exact solutions for vacuum Bianchi models in JBD [9, 10, 12]. 

In the reparametrized formulation of the equations for the anisotropic homo- 
geneous metric of Bianchi type VIIq, exact solutions can be obtained for the case 
of a Bianchi- type VIIq; the behaviour of such solutions depends on the sign of the 
quantity 

A = -4(^ + 1/4), (12) 

where, as can be readily shown, 

B^^'Hl + -Hi-—H[, (13) 

is a constant, i.e. it is a first integral of the system (10) and (11) expressed in terms 
of the Hubble expansion rates. It can be shown that, out of the possible solutions 
of equations (10), the only physically plausible is the one corresponding to the case 
A < 0; the cases A > or A — can be shown to led to negative or complex scale 
factors. For some details see [9]. The only physically admissible solution can then 
be explicitly written as 

ai($)=(^^-^J ($cosh [-V(S + l/4)ln(/$2)]) ^ (^4) 
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where cq and / are positive integration constants. The other two scale factors can 
be easily obtained from ai($) and from equations (6) and (10), they are 



= coai($), a3($) = 2-^/^coai{^). 



(15) 



These results show that the three scale factors are proportional to each other. This 
means that the Bianchi VIIq model, despite what we could have anticipated, shows 
an isotropic evolution; it also implies that the shear a, vorticity f2 and acceleration 
a of the reference congruence, all vanish. The vacuum Bianchi- VIIq JBD universe 
thus basically behaves as Friedmann- Robertson- Walker (FRW) space-time [14]. See 
figure 1. The local volume on the surface of homogeneity is 



The constriction equation implies the following relationship in our case 



from here we notice that to have meaningful solutions the coupling parameter has 
to be restricted to uj > —3/2 = Um or, B > —1/4 ; oOrn is just the minimum value uj 
could have for the JBD theory to make sense (Ruban and Finkelstein 1975). Notice 
that the specific value of uj is enough to determine the evolution of the Hubble 
expansion rates through equation (13), whose solution is (assuming B > 0) 



where Hq = iJi($o)- As the 3 scale factors are proportional to each other, it then 
follows that the three Hubble rates are always the same: Hi = H2 = H3 — another 
manifestation of the equivalence to a FRW isotropic universe. From (18) it can be 
seen that the Hubble rates, after reaching a minimum value, vanish asymptotically 
as ~ — log($)/$ when $ ^ 00. 

Using equations (9) and (16) we can obtain the dependence of $ on t, as follows 



V = 




(16) 




(18) 
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(AB + 1)3/^ 
V2co 




$cosh[VB + l/4 ln(/$2)] 



3/2 



(19) 



though we can obtain t as a function of as in (19), we cannot invert it to obtain 
exphcitly $ as a function of t. Nevertheless, it is easy to realize the enormous change 
that occurs in $ over a very small span of t values. Behaviour of this sort has been 
shown to have relevance in explaining, for a Bianchi IX universe, the vanishing of 
the maximum Lyapunov exponent calculated in the intrinsic time $ in spite of it 
being positive in the synchronous time t [13]. This also shows that, asymptotically, 
$ grows without bound whereas t approaches a certain finite value te- As follows 
from (19), this particular value of the cosmological time can be evaluated as 



where we defined a = 3(^-B + 1/4 - 1/2), b = A^B + l/A. Under the conditions of 
the Bianchi model we are working with, expression (20) has always a i?-dependent 
finite limit — which can be expressed in terms of the so-called Barnes's extended 
hypergeometric function [15]. See figure 2. This means, of course, that in order to 
analyse the fate of the Bianchi- VIIq universe for cosmological times greater than 
a different choice of coordinates is needed. 

We can now state the consequences of (14) for the behaviour of the universe. 
An analysis of the expression (14) shows that, if a; > (hence, only if B > 0), 
the universe begins with a singularity since lim$^o o^i = 0, see equation (22), and 
then, as the intrinsic time unfolds, the universe expands as ~ where a = 
3^B + 1/4 - 3/2; the universe then reaches a maximum volume Vmax; that can 
be easily calculated from (14) and (16), and then shrinks until it reaches again a 
singularity. In the other case (B or u < 0), as it can be ascertained from (14), (16), 
and (22), the fate of the universe depends of whether uj is larger than Uc = —4/3 
or not. If it is larger than Uc (though still negative), it is possible to show that 
lim<j,_jO (^{t) — which coupled with the expression for the Ricci scalar (22) means 
that the universe starts from a non singular, spatially infinite and fiat {R = 0) state, 
and from there it isotropically collapses to a singularity; that is, lim^^co o^(^) — 0) 
lim$_>oo R — oo. 




(20) 
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A universe is singular if the value of the Ricci scalar R — g°''^Rab along a certain 
geodesic congruence blows up, i.e. if i? ^ ±00, whereas the associated affine 
parameter s tends to a finite value [16]. If this happens, then we say that there 
exists a curvature singularity and thus that we are dealing with a singular spacetime. 
For our Bianchi VIIq model, we have chosen as the congruence the world lines of 
test observers (time-like geodesies) whose affine parameter is the synchronous time 
t. The scalar curvature can be readily shown to be [9] 



using only equations (5), (10) and (11) we can rewrite R for any vacuum Bianchi- VII 
model in terms of the coupling parameter uj and the scalar field 0, as follows 



It is important to notice that expressions (21) and (22) are valid for all the Bianchi- 
type VII models, not just for the specific VIIq case [9]. In order than i? ^ 00 in (21), 
all that is needed is that at least one the scale factor ai, or just the (/)-field, vanish at 
a finite value of the synchronous time t; that is, that they vanish for whatever value 
of the intrinsic time $. Notice also that we can regard the evolution of the Bianchi- 
type VII universe — equations (3) and (4) — as "driven" by the curvature R. An 
analysis of the solution (14) together with (22) (basically expanding in series the 
expressions for ai($),) shows that, the possible singularities that may arise depend 
on the value of according to the following Table 1. 

The role of the curvature scalar in governing the universe can be qualitatively 
described as follows (Table 1): a) If u; > 0, a strong curvature prevents expansion; 
it is only when curvature is small that the full extent of expansion is reached 
but, as soon as the curvature is large in magnitude again, contraction sets in. If 
> a; > oJc, the strong curvature of the initial singular but infinite in extent state 
initiate the contraction until R is very small, after that the universe collapses as 
$ ^ 00 reaching a localized singularity. If > u; > —3/2, the universe begins in a 




(21) 




(22) 
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non-singular flat infinitely extended state and from there it collapses to a localized 
singularity as $ ^ oo. 

Let us conclude by saying that the supposed homogeneous and anisotropic 
Bianchi VII model in fact shows an isotropic expansion in the case in which h = 0. 
On the other hand, we show that the dynamics of the early stages of the expansion 
in this JBD model mainly depends on just one of the scale factors (that, here, we 
choose as ai) and that much of the universe behaviour depends on the value of uj. 
We have also obtained the dependence of the three scale factors ai on the intrinsic 
time $. As we have concluded that the three scale factors are proportional to each 
other, the expansion is necessarily isotropic and not anisotropic, as it is usually 
assumed to be in this model. In fact, we have shown that a Bianchi- VIIq JBD 
vacuum universe is basically equivalent to a FRW-spacetime with a — 0, ft — and 
a = 0. 

We think an important point of our analysis is that shows that, even starting 
with supposedly anisotropic models, the inclusion of a scalar field can drastically 
isotropize the behaviour offering the possibility of coordinating the model with 
the observed isotropic properties of the actual Universe. Moreover, our analysis 
shows how we can change the dynamics by just changing u, a type of behaviour 
which gives support to the extended infiation idea of using a dynamical coupling 
parameter [4, 13]. It should be clear also that from the coordinates used in this work, 
we cannot obtain information about what happens with the model for synchronous 
times greater than te, (for example, ^what happens as i — > oo?). From this point 
of view, the analysis performed refers only to the early epochs of the universe and 
thus it just describes its early fate. 
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Figure Captions 

Figure 1 

The behaviour of the three scale factors ai, 02 and 03 in a Bianchi-type VIIq 
universe against the intrinsic time is shown; these are plots of equations (14) and 
(15). Notice that the 3 scale factors are proportional to each other; the behaviour 
shown is generic in each interval of u values. 

a) The case a; > 0, the especific values used for the parameters are B — 33.333, 
Co = 2, / = 1. 

b) The case > a; > Uc, the values used for the parameters are B — —0.067, cq = 2, 

c) The case Uc > uj > —3/2, the values used for the parameters are B — —0.225, 
Co = 2, / = 1. 

Figure 2 

The maximum cosmological time tg (equation (19)) for which the analysis gives 
information against the value of B{— u>/6). 

Figure 3 

The absolute value of the scalar curvature \R\ in the Bianchi-type VIIq universe 
is shown in the range < |i?| < 00 versus $ in the range < $ < 00. What we 
really graph here is arctan(|i?|) against arctan($). Small ^-values correspond to 
small t-values but large values of $ correspond to t = te — 3.624, as can be seen in 
figure 2. Giving the way the graph was compactified, you should interpret as 0, 
but 7r/2 as 00. 

a) The case a; > 0, the especific values used for the parameters are B — 33.333, 
Cq = 2, / = 1. The universe goes from a point singularity and end in another. 

b) The case > a; > lOc, the values used for the parameters are B — —0.067, cq = 2, 
/ = 

c) The case c^c > > ujm, the values used for the parameters are B — —0.225, 
CO = 2, / = 1. 
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Table caption 

Table 1. 

The behaviour of the Bianchi VIIo model for the different uj values is summarized 
in this table. R is the scalar curvature. The critical value of the coupling parameter 
is ujc — —4/3. For uj < uj^ = —3/2 the JBD theory loses sense. 
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Coupling Parameter 


Initial Size 


Final Size 


Initial R 


Final R 


u; > 








— oo 


— CO 




finite 











0>u;> -4/3 


oo 





oo 


oo 


uj = -4/3 


oo 








oo 


-4/3 > a; > -3/2 


oo 








oo 
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